Collimation of energy in medium-modified QCD jets 
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Abstract: The collimation of energy inside medium-modified jets is investigated in the leading log- 
arithmic approximation of QCD. The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equations are slightly modified by introducing splitting functions enhanced in the infrared sector. As 
compared to elementary collisions in the vacuum, the angular distribution of the jet energy is found to 
broaden in QCD media. 
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1 Introduction 



The jet quenching has been established as the phenomenon of strong high-pr suppression of partons 
produced in heavy ion collisions 0~H3]]. Experimentally, it has been confirmed that highly virtual partons 
produced in such reactions suffer an energy degradation prior to their hadronization in the vacuum. The 
high-pr partons energy is radiated away as gluon bremsstrahlung, thus increasing the amount of such 
gluons after several interactions of the leading parton with the scattering centers in the quark gluon 
plasma (QGP) occur 10). Jet production in hadronic collisions is a paradigmatic hard QCD process. An 
elastic (2 — > 2) or inelastic (2 — > 2 + X) scattering of two partons from each one of the colliding hadrons 
(or nuclei) results in the production of two or more partons in the final-state. At high pt , the outgoing 
partons have a large virtuality Q which they reduce by subsequently radiating gluons and/or splitting into 
quark-antiquark pairs. Such a parton branching evolution is governed by the QCD branching probabilities 
given by the DGLAP equations [5] down to virtualities 0(1 GeV 2 ). 

In this paper we study the energy collimation of jets produced in heavy ions collisions from a toy QCD- 
inspired model introduced by Borghini and Wiedemann JQ. In this model, the DGLAP splitting functions 
are enhanced in the infrared region in order to mimic medium-induced soft gluon radiation. More pre- 
cisely, the 1/x dependence of the branching probabilities q — ^ qg, g — > gg is enhanced by introducing 
the factor N s = 1 + f me d in P a b(x) = (1 + fmed)/% + 0(1). In order to quantify the collimation of 
energy, we consider a jet of half opening angle Go initiated by a parton A, followed by the production of 
a subjet of half opening angle Q (with G < Go) initiated by a parton B with large energy fraction x ~ 1, 
which contains the bulk of the jet energy registered by the calorimeter. Experimentally, this corresponds 
to the calorimetric measurement of the energy flux deposited within a certain solid angle Q. The first 
version of this quantification was given in QUI in the vacuum and it was found to scale like, 



where 7,4(2;) (A = q, q, g) is a certain function derived from the DGLAP evolution equations at large 
x and A = Aqcd is the mass scale of QCD. In this context, we will introduce the N s dependence in 
7(x) — > -y(x,N s ) in order to compare our results in the medium (j(x,N s )) with previous results in 
the vacuum ('j(x) = j(x, 1)). The same logic was applied to the study of the collimation of average 
multiplicities in |9 ], where jets were found to broaden as compared to those in the vacuum after inserting 
the N s -dependence in the evolution. 

The same model has been largely used, i.e. for the estimation of other observables like inclusive px- 
distributions ifTOTl and fragmentation functions iTTTTl . in all cases results look promising and appealing for 
future comparison with experimental data. The large suppression of these distributions II 1 Oil 1 1 H at large 
x are good examples of the large energy suppression suffered by high-py partons in the dense QGP. 

2 Theoretical framework 

Let us start from the production of one gluon or quark (A = g, q, q) initiated jet with energy E, followed 
by its fragmentation into another parton B at angle Q such that, Qo/E < G < Go, where Go is the half 
opening angle of the jet and Qq, the minimal transverse momentum of the emitted partons before the 
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hadronization into hadrons occurs. The collimation of the energy is characterized by the energy fraction 
x ~ 1 where the bulk of the jet energy inside the given cone < ©o <C 1 is deposited. The probability 
for the energy fraction x to be deposited in a cone of aperture G is related to the inclusive spectrum of 
partons through the formula [8 ], 



D A (x,EQ ,EQ) = D^(x,EQ ,EQ), 

B=g,q 



(1) 



where the nature of partons B is not identified. In the above relation dH) two scales EQq and £0 have 
been written in the argument of the inclusive spectrum so as to account for the given process in the energy 
range EQ < Q < EQq; this notation will be shortened below. In order to quantify the collimation of 
energy, it should be considered that the deposited fraction of the jet energy is large, that is for x — >■ 1. 
The DGLAP evolution equations in the vacuum takes the simple form 0, 
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;D(x,Q 2 



a s (Q 2 
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dz 



P(z)D(-,Q 



(2) 



where P is the Hamiltonian matrix of splitting functions characterizing the parton splitting probabilities. 
The splitting functions are accurately obtained from the FO approach in perturbation theory for small 
coupling constant q s <1. The coupling constant is written in the form, 



a s (Q 2 



7T 



NcPoHQV^qcd) 



where (3q = jj^(^-N c — is the first coefficient of the beta function and Aqcd ~ 250 MeV. In 
Mellin space V(j, Q 2 ) is obtained from the transformation 



V(j,Q 2 )= [ dxx j - l D{x,Q 2 
Jo 



such that the convolution (O simply becomes 

d 



dlnQ 2 

or more explicitly rewritten in the matrix form at LO 



nj,Q z ) = r(j)v(j,Q- 
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V qq {j) V qg (j) 
Vggti), 



VMS) 

V Vg(j,0 J 



(3) 



where T> qNS and V qs stand respectively for the flavour non-singlet (or valence) and flavour-singlet quark 
distributions, and Vik(j) is the Mellin transform of the LO splitting functions. The variable Q 



e(Q 2 



In In 



Q 1 



A 2 



dlnQ 2 



a s d 



is introduced for the sake of simplicity. In order to account for the medium-induced gluon radiation in 
heavy-ion collisions, we make use of the QCD-inspired model proposed in (5), which allows for a simple 
computation of the equations at large x. In this model the infrared parts of the splitting functions are 
enhanced by the factor N s = 1 + f me d, where f me d > accounts for medium-induced gluon radiation. 
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Realistic values of f me( i (= 0.6,0.8) are extracted from fits of the nuclear modification factor -R^A^to 
the RHIC data [2]; for the LHC we will take f me d = 1. The medium-modified splitting functions read, 



Pgg{z)=W e 



P qg (z) = 2C f 
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, P gq (z) = 2T R [z 2 + (l-z) 2 } 



2N S 
1 - z 
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(4a) 
(4b) 



with the [...]+ prescription defined as L 1 dx[F(x)]+g(x) = J dxF(x)[g(x) — g(l)]. Performing the 
Mellin transform of Eq. (I4al4bb gives iTTTIl 



P 99 U) 
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(5a) 
(5b) 
(5c) 
(5d) 



where if? (J) = InT(j) is the digamma function with ~ 0.5772 the Euler constant. It can be easily 
checked that Eq. ((5]) reduces to the ordinary splitting functions given in J5HH after setting N s = 1. In 
the large z ~ 1 limit, or equivalent large j > 1 we are interested in, the expressions of the anomalous 
dimensions (|5al > and (f5d]) can be re-expressed in the form 



V qq (j) ^4C F N S [- In j + 



4iV s 



Je), P gg (j)~W c N s -lnj + 



iV, 



IE 



(6) 



where the asymptotic behavior if)(j + 1) sj lnj has been replaced for j>l. The Mellin transform can 
be inverted, with the inversion given by 



D(x,A£) 



2iri 



djx-W(j,A£), 



(7) 



c 



where the contour C in the complex plane is parallel to the imaginary axis and lies to the right of all 
singularities of A£). In ^} 
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Integrating (0) after inserting © leads to the valence distribution at large x ~ 1 Q, 



Di(x,A0 ~ (1 - x )-i+4C^ s a ? 



exp[4C A jV s ( g | rr - 7g )Ag] 
T(4CUiV s A£) 



(8) 



where /?o = 3/4 for n/ = 3 was replaced from © and the N s dependence is new. The behavior of 
the valence distribution ([8]) when 6 changes can be studied by taking its derivative with respect to £. 



3 Raa corresponds to the ratio of medium-modified and unmodified single inclusive hadron spectra. 
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In (SI, the valence distribution ([8]) was shown to present a maximum at some angle for fixed values of 
x where the bulk of the jet energy is concentrated. Taking d^D^ = results in, 
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^(4CaN s A£) = ln(l - x) + 



AN, 



IE, 



(9) 



where here again ij)(x) is the digamma function. For A£ — >• 0, or — )■ ©o, almost the whole energy will 
be deposited inside the cone ©o while for decreasing down to © > E/A the emission outside the cone 
increases. The equation (© with the N s dependence is the main result of this paper. It should be inverted 
numerically in order to provide the behavior of the ratio @/@o as a function of the energy of the leading 
parton and the nuclear factor N s for fixed values of x ~ 1. Symbolically, the inversion can be written in 
the simple form, 
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(10) 



such that, the collimation is stronger for higher energies in both vacuum and medium jets. As observed 
in (PTOl . the N s -dependence contained in the denominator of the exponent should soften the collimation 
of jets produced in heavy ion collisions. Furthermore, the function ja( x , N s ) appearing in the exponent 
provides the slope of the collimation as a function of x and N s . 



3 Phenomenology 

We display the collimation of energy through the ratio from the numerical inversion of (©. Let us first 
consider a jet produced in a heavy ion collision at energy scales greater than 100 GeV at the LHC and let 
us set ©o ~ 1 for simplicity. For the gluon and quark jets, one sets Ca = N c = 3 and Ca = Cp = 4/3 
respectively in (©. In figured] and |2j we display the collimation of energy for two fixed values of energy 
fractions x = 0.5 and x = 0.9 which are in agreement with the large x (x ~ 1) approximation applied 
in this frame. As expected from Q and observed in figures Q] and |2j the curves of the collimation in 
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Figure 1: Collimation of energy inside a gluon jet for x = 0.5 (left) and x = 0.9 (right). 



the medium (f me d = 0.6, 0.8, 1) are shifted towards higher values of the subjet opening angle as 
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compared with the curve for in the vacuum. Moreover, the slope of the collimation decreases as f me d 
increases. For fixed @o, E and x, the same amount of energy contained inside a subjet of opening angle 
61 in the vacuum is distributed over a broader angular aperture 02, such that ©2 > ©1. It follows 




Figure 2: Collimation of energy inside a quark jet for x = 0.5 (left) and x = 0.9 (right). 

that gluon jets are less collimated than quark jets in the vacuum [8] and not surprisingly, also in the 
medium. For x = 0.5 jets are more collimated than for higher values of the energy fraction. Finally, it 
is straightforward to check that the values of 7,4 (ar, N s = 1) for x = 0.5, 0.9 are in agreement with the 
expected values 7 9 (0.5, 1) » 0.54, 7 9 (0.9, 1) « 0.30, 7 9 (0.5, 1) « 0.83 and 7 9 (0.9, 1) « 0.55 in the 
vacuum 0[H. This analysis could be performed in future measurements by the ALICE collaboration as 
a probe for the QGP produced in PbPb-collisions. 

4 Summary 

In this paper we have studied the collimation of energy inside jets produced in heavy ion collisions. From 
the production of one jet of opening angle ©o initiated by a parton A, we considered the production of 
a subjet (0 < ©o) initiated by a parton B, where the definite fraction x ~ 1 of the jet energy is 
deposited. We made use of the DGLAP evolution equations at large x with medium-modified splitting 
functions, which accounts for medium-induced soft gluon radiation. Our study was extended to the 
distributions D^(x, Q 2 ) at x ~ 1 so as to provide their dependence on the nuclear parameter N s and 
therefore, it allowed to determine the behavior of the collimation of the jet energy as a function of the 
same parameter. 

As physically expected, jets in the vacuum are more collimated than medium-modified jets. Hence there 
is an evidence for the broadening of jets as gluon radiation is enhanced. Quark jets are more collimated 
than gluon jets in both the vacuum and the medium. Therefore, the collimation of energy treated in this 
paper and the collimation of average multiplicity treated in [9 ] are both good candidates to investigate 
the jet quenching in heavy ion collisions. Finally, more efforts are required towards the construction of a 
more realistic QCD picture where the medium parameters and medium-induced soft gluon dynamics lfi"2l 
could be taken into account. 
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